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Abstract 

The purpose of this paper is to present an algorithm for evaluating Hankel 
transform of the null and the first kind. The result is the exact analytical 
representation as the series of the Bessel and Struve functions multiplied to 
the wavelet coefficients of the input function. Numerical evaluation of the 
test function with known analytical Hankel transform illustrates proposed 
algorithm. 

The Hankel transform is a very useful instrument in a wide range of physical 
problems which have an axial symmetry pp . The influence of the Laplasian on a 
function in a cylindrical coordinates is equal to the product of the parameter of the 
transformation squared and the transform of the function: 
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The Hankel transform of the null (n = 0) and the first (n = 1) kind are repre- 
sented as 

oc 

F n(p) = J f(r)J n (pr)rdr, 

So (2) 
fn(p) = J F(p)J n (pr)pdp. 



Besides that integrals like (2) are connected with the problems of geophysics 
and cosmology, for example 

But practical calculation of direct and inverse Hankel transform is connected 
with two problems. The first problem is based on the fact that not every transform 
in the real physical situation has analytical expression for result of inverse Hankel 
transform. The second one is determination of functions as a set of their values for 
numerical calculations. Large bibliography on those issues can be found in The 
classical trapezoidal rule, Cotes rule and other rules connected with replacement 
of integrand by sequence of polynoms have high accuracy if integrand is a smooth 
function. But f(r)J n (pr)r (or F p (p)J n (pr)p) is a quick oscillating function if r 
(or p) is large. There are two general methods of the effective calculation in this 
area. The first is the Fast Hankel transform The specification of that method 
is transforming the function to the logarithmical space and fast Fourier transform 
in that space. This method needs a smoothing of the function in log-space. The 
second method is based on the separation of the integrand into product of slowly 



1 



varying component and a rapidly oscillating Bessel function [7|- But it needs the 
smoothness of the slow component for its approximation by low-order polynoms. 

The goal of this article is to apply wavelet transform with Haar bases to . 

The both direct and inverse transforms J2J) are symmetric. Let us consider only 
one of them, for example, direct transform. Let's denote f(r)r as g(r). Then Hankel 
transform is 



FoAp) = J g(r)Jo,i(pr)dr (3) 
o 

The expansion g(r) 6 L 2 (R) into wavelet series with the Haar bases is [S] : 

oo oo oo 

9 (r) = c ° k ^ k M + Y Y d i k ^jk ( r ) . ( 4 ) 

fc=0 j=0 k=0 
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(r) = V H (r - k),i> jk = V 2 ^ H (Vr - fc), 



1, te (0,1/2) 
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After substitution (J2J into © one has 

oo °? oo oo °? 

cok / yfc(r)Jo,i(pr)(ir + y^y^rfjfc / ip jk (x)Jo,i [pr)dr. 
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Making use integrals from 2 we have as a result 
F (p) = E c ofe [(fc + l)Jb (p(* + 1)) - kJ (pk) 

keZ 

f [(fc + 1)D (p(k + 1)) - kD (pk)} + 



E E djfc [2(* + 5)^0 (p(* + I)) - (k + 1) Jo (p(fc + 1)) - fc Jo (pk) - 
j=o kez 

§ [2(fc + \)D (p(k + I)) - (fc + l)D (p(k + 1)) - kD (pk)] 

Fi (p) = i ( E cofe [Jo (pfc) - Jo (p(fc + 1))] + 

E E rfifc [2Jo (p(k + |)2-J) - J (p(fc + 1)2- J ) - Jo (pAriT*)] 
j=ofcez 



(5) 
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where -D(£) = -Ho(£) Ji(£) ~ ^l(S) Jo(£) an d -ffo.i is Struve function of the null 
and the first kind. 

The most sufficient result is that equations iJSJ and © are exact. They can be 
used in any analytical expressions. Especially it's useful for Hankel transform of the 
first kind because © contains only a combination of Bessel functions, and one can 
use such their properties as orthogonality etc. The coefficients cok means average 
value of g(r) at the range [fc, fc + 1] : 

k+l 

c k = J g{r)dr 

k 
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Figure 1: Original function (left) and Hankcl transform (right) 
The detail coefficients are 

2- 3 {k+i) 2~'(fc+l) 

= 2-' 1 \ J g(r)dr - J g(r)dr 

The formulas an d © allow us to get full analytical solution if integrals above 
have close form solution. In the opposite case the solution must be numerical but 
this method provides an effective algorithm for that. It's obvious that djk decrease 
very quickly if g(r) is a smooth function. One can practically use djk > e, where 
e is small. But if g(r) has steps, sharp vertices or discontinues then the detail 
coefficients concentrate around these points and one can appropriate they are equal 
to the zero in other areas. 

Let us consider for example a function with known analytical Hankel transform 

e^^rJxiprydr = JL e -P , /4«\ (7) 
4a 4 
o 

The approximation and detail coefficients may be calculated analytically in closed 
form 
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Thus © with the coefficients (JSJ is the exact representation of the Hankel trans- 
form. Let us consider the approximate solution. Suppose the function is known l(7|l 
only in the segment [0, h]. Then there is the series instead of 10} : 

J 2 J -1 

g(r) = Co<p (r) + ^ d j^jk( r )- (9) 

If J — > oo then is exact for this truncated function. But practically we only 
use several first levels. For example we can see the original function (the replacement 
r to x = r/h is used) (J7J) and the transform at the Fig.l. One can see that exact 
transform (solid line) and the transform at level J — 3 (dotted line) coincide at 
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this figure. The absolute errors between the exact transform and the approximate 
transform at the levels J = 2 (solid line), J — 3 (dashed line), and J = 4 (dotted 
line) are represented at the Fig. 2 (left). It's oblivious that the error is small in 
comparison with the values of the F\{p). The absolute error at the level J = 3 in 
a wide range of p is plotted in the right side of the Fig. 2. One can see that this 
error has quasi-periodic oscillations because the function is truncated. But they 
decrease with the growth of p (and J) when oscillations in classical Fast Hankel 
Transform 4 increase. 
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